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Abstract. - Non-equilibrium structure formation and conversion of spinning to translational 
motion of magnetic colloids driven by an external rotating magnetic field in microchannels is 
studied by particle-based mesoscale hydrodynamics simulations. For straight channels, laning is 
found. In ring channels, the channel curvature breaks symmetry and leads to a net fluid transport 
around the annulus with the same rotational direction as the colloidal spinning direction. The 
dependence of the translational velocity on channel width, ring radius, colloid concentration, and 
thermal motion is predicted. 
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Externally actuated and self-propelled micro- and nano- 
rotators show an intriguing variety of non-equilibrium 
structure formation and dynamics. Examples of such sys- 
tems include super-paramagnetic colloidal particles in a 
rotating magnetic field [IHl] , dipolar colloids in a rotating 
electric field [5] , colloidal dimers rotated by laser tweezers 
[5] and biological swimmers such as Volvox algae [7] . The 
collective behavior of rotator suspensions is governed by 
hydrodynamic interactions. While a single rotator can- 
not show any directed translational motion, in ensembles 
of several rotators their spinning motion leads to an en- 
hanced translational diffusion [8^. Other interesting phe- 
nomena in suspensions of rotators are stable bound states 
of spinning Volvox algae [7], spontaneous pattern forma- 
tion of spinning magnetic disks at liquid-air interfaces in 
the form of rotating hexagonal crystal structures [T] , and 
spinning colloids, placed asymmetrically in a microfluidic 
channel, which act as micropumps [3J|3J[B]. 

A key prerequisite for the conversion of rotational into 
translational motion is symmetry breaking in confinement. 
For a single rotating cylinder close to a planar wall (and 
oriented parallel to the wall), the resulting hydrodynamic 
forces vanish since the viscous stress is balanced by the dy- 
namic pressure field [S]; thus, a single wall is not sufficient 
to generate translational motion of a cylinder (in contrast 
to a sphere, which experiences a force parallel to the wall 
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in the direction expected for a rolling motion [10]). How- 
ever, in the presence of a second parallel wall, more fluid is 
pumped through the wide gap between cylinder and wall 
than through the narrow gap, causing a net flow in the 
channel and a reaction force parallel to the wall on the 
cylinder — unless it is centered between the walls [TT1[T2] . 
Pumps based on this principle, with a typical channel di- 
ameter in the centimeter range and a fixed position of 
the rotating cylinder, have been studied both experimen- 
tally and numerically [TTJ[T2]. The goal of microfluidics 
is to work with much smaller length scales, on the order 
of micrometers or smaller. Indeed, micron-sized pumps 
employing rotating colloids have been contructed recently 

I3I1E]. 

With further miniaturization of microfluidic systems, 
thermal fluctuations become increasingly important. Fur- 
thermore, a sufficiently large concentration of particles is 
required to achieve reasonable flow velocities. Thus, we 
investigate a system of spinning colloidal particles, which 
are free to move in a channel, are small enough to dis- 
play thermal Brownian motion, and have a finite concen- 
tration. Pairs of co-rotating cylinders in an unbounded 
viscous fluid are known to mutually exert forces on each 
other that are perpendicular to their connecting line |13j . 
which implies a circling motion. At higher concentrations, 
coordinated motion occurs [5]. 

In a channel, many-body effects and the presence of the 
walls modify the interactions, and the situation becomes 
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Fig. 1: (Color online) Simulation snapshots of rotating colloids 
in a ring channel (upper part) with median radius R = 2.4(7 
and width D = 2. Oct, where cr is the colloid diameter, and of 
a section of a linear channel (lower part) with D = 2.5(7. The 
rotational and translational motion of the colloids is illustrated 
by the arrows. For movies, see also footnote 1. 



much more complex. Symmetry implies that there cannot 
be any net fluid or colloid transport in a straight channel. 
However, we find that the particle distribution across the 
channel causes local directed motion and leads to lane for- 
mation. In ring channels, the additional symmetry break- 
ing leads to a net translational motion along the channel, 
as illustrated in Fig. [TJ Simulation animations of spin- 
ning colloids in a linear channel (linear_channel.avi) 
in a circular channel (ring_channel.avi) are available in 
the online supporting material. We address the funda- 
mental question how symmetry breaking due to channel 
curvature generates net transport in the presence of ther- 
mal fluctuations. For simplicity, we restrict ourselves to 
two-dimensional systems. 

The fluid is modeled by multi-particle collision dynam- 
ics (MFC), a mesoscale hydrodynamics simulation tech- 
nique that naturally includes thermal fluctuations jl4H16j . 




^ The movie ring_chaiinel.avi shows a simulation animation of 
2D colloids in a ring channel of median radius R = 2Act and width 
D = 2.0(7. The resulting fluid velocity field is displayed by arrows. 
The movie linear .channel . avi shows counter-clockwise rotating 2D 
colloids in a linear channel of width D = 2.5(7. In the displayed 
system, the colloid area fraction is <& = 0.38. In both movies, a 
constant external torque Lcxt = ISOfcaT is applied. 



Fig. 2: (Color online) Colloid density distribution across a ring 
channel of median radius R — 2.4a and width D = 2.5(7, with 
area fraction $ = 0.25. Data are shown for torques Lcxt,o (full 
red line) and I/ext,o/3 (dashed green line), as well as for colloids 
at thermal equilibrium (dotted blue line). Inset: Same for a 
straight channel with torques Loxt,o (full red line) and Lcxt,o/9 
(dashed green line). 



Here, the fluid is represented by point particles; their dy- 
namics proceeds in two alternating steps. In the streaming 
step, the fluid particles move ballistically. Subsequently, 
they are sorted into the cells of a square lattice with lat- 
tice constant a. New relative velocities (with respect to 
the centcr-of-mass velocity of each cell) are assigned to 
all fluid particles in a collision cell, which mimics the si- 
multaneous interaction of these particles. The algorithm 
is constructed such that linear momentum is conserved 
in each cell, in order to obtain correct hydrodynamic be- 
havior. For a study of hydrodynamics of rotating colloids 
with fixed torque, it is essential to employ a variant of 
MFC which also conserves angular momentum [TTHH]. 

The colloids are modeled as discs of diameter a with no- 
slip boundary conditions. For the colloid-colloid and the 
colloid-wall interactions, a purely repulsive (shifted trun- 
cated) Lennard- Jones (LJ) potential is added. Its range 
Slj ~ O.ltT is chosen to be small compared to the col- 
loid radius, but large enough to ensure that there is at 
least one collision box of fluid between surfaces. This 
guarantees that hydrodynamic interactions are captured 
quantitatively even at the smallest colloid distances [18]. 
We apply a constant external torque Lext to all colloids, 
causing the colloids to spin counter-clockwise about their 
axes with an average angular velocity ftoo — Lcxt/iwo-'^) 
in the bulk, where rj is the viscosity of the fluid. In 
a channel, the colloids develop a characteristic transla- 
tional velocity vq. The Peclet number for the rotating col- 
loids is then defined by Pe = avo/Dt, where Dt is the 
(geometry-dependent) translational diffusion coefficient. 
The Peclet number depends both on the channel geometry 
and the external torque. For the largest torque consid- 
ered, Lext = ISOfcsT, the Peclet number varies between 
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Fig. 3: (Color online) (a) Colloid velocity profiles for various 
colloid area fractions $ for a straight channel of width D — 3a. 
(b) The corresponding fluid velocity profiles. 



Pc = 0(10) (narrowest channel) and Pc 0(10^) (widest 
channel). For the widest channel, the Peclet number varies 
linearly with the colloid torque, with Pe = 0(10) for the 
smallest torque £cxt,o/9 = IG.TfcsT. If not mentioned 
otherwise, the torque Lext,o is employed. 

First, we investigate straight channels that are wide 
enough to allow colloids to pass each other, with widths 
D = 2.5(7 and D = 3a. We observe that the spinning 
colloids move in opposite directions at opposing walls (see 
movie linear .channel . avi), and with increasing concen- 
tration approach the walls very closely, see Fig. [21 with a 
distance only slightly larger than a/2 + A compar- 
ison with simulations for smaller spinning velocities, see 
Fig. [5J shows that the peaks in the density profile broaden 
with decreasing Lext- A weak density increase near the 
walls due to the finite colloid concentration remains even 
in thermal equilibrium. For wider channels, the surface 
excess is less pronounced. 

On a qualitative level, this behavior can be understood 
as follows. When two spinning colloids pass each other in 
the channel, the hydrodynamic forces mutually exerted on 
each other enhance this directed motion. For symmetry 
reasons, colloids spinning in the center of the channel can- 
not propel themselves forward; however, they also cause a 
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Fig. 4: (Color online) Average colloid velocity per half-channel 
as a function of the colloid area fraction $ for straight channels 
of width D = 2.5(7 (□, ■, red) and D = 3a (o, •, blue). Data 
are shown for three different torques, Lcxt.o (B, •), icxt,o/3 (□, 
o), and icxt,o/9 (inset). 



local fluid motion, and hence push nearby colloids towards 
one of the walls, while simultaneously being pushed to the 
opposite wall. 

Colloid velocity profiles across the channel are shown in 
Fig. 13^ for Z? = 3(7 and various colloid densities, together 
with the corresponding fluid velocity profiles in Fig. 
The fluid velocity profile is a result of the combination 
of colloid motion pushing the fluid column ahead of it 
forward, the surface velocity due to colloid rotation, the 
no-slip boundary conditions on the walls, and the veloc- 
ity field generated by the colloids near the opposite wall. 
Moreover, since a constant torque is exerted on the col- 
loids, they spin faster at larger distances from the wall 
and at smaller (local) colloid densities. Note that the di- 
rected motion near the walls does not imply that individ- 
ual colloids always move in the same direction; instead, 
they change lanes randomly due to thermal fluctuations 
and hydrodynamic interactions (which promote a circling 
motion of near-by colloids). 

Although the fluid and colloid velocities averaged over 
the whole channel vanish due to symmetry reasons, there 
is a directed motion in both half-channels. To quantify 
this effect, we calculate the average colloid velocity per 
half-channel, as shown in Fig. |4l In the limit of vanishing 
colloid concentration, the only forces acting on the col- 
loid originate from the hydrodynamic interactions with the 
walls. In contrast to Fig. [51 the density distribution across 
the channel for a single colloid is a uniform distribution for 
small Pe; because there is no propulsion near the channel 
center, but a propulsion in opposite directions near oppos- 
ing walls, the velocity per half-channel approaches a finite 
but small value in this limit (see Fig. [4]). With increasing 
colloid density, pairwise interactions between colloids be- 
gin to contribute. As mentioned above, when two colloids 
come close to each other, they circle around each other. 
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Fig. 5: (Color online) Average tangential velocity i;tan of all 
colloids as a function of the median annulus radius R for curved 
channels with widths D = 1.5(7 (•, red), 2. Oct (■, green), 2.5a 
(a, blue) and 3. Oct (□, magenta). The torque on each colloid 
is Z/cxt,o. 



Fig. 6: (Color online) Schematic illustration of colloids spinning 
close to the inner and outer confining cylinders of a ring chan- 
nel. Fluid is pumped towards or away from the walls (curved 
black arrows), giving rise to dynamic pressure. Excess pres- 
sure is indicated in yellow (light gray), low pressure in red 
(dark gray). The small straight white arrows depict the shear 
flow in the gaps. 



0.1 



a, 



0.01 



which implies an increased and reduced density near the 
walls and in the center, respectively. Because the propul- 
sion is most efficient near the walls, the colloid velocity 
in the half-channels increases when the colloids pass each 
other, and the fluid volume each colloid has to move de- 
creases as 1/$, the half-channel velocity increases with $ 
at low densities. 

At very high colloid concentrations, the average colloid 
velocity also becomes small because lubrication forces be- 
tween neighboring colloids and between colloids and the 
walls slow down the rotational motion. Hence, the half- 
channel velocity exhibits a maximum velocity at a finite 
colloid area fraction $ ~ 0.15, see Fig. Moreover, the 
velocity is higher in the broader channel, where the col- 
loids can pass each other more easily. 

Fig. |4]also demonstrates the effect of a variation of the 
colloidal spinning velocity by a change of the external 
torque. Wc observe that Whaif/(o'i^oo) is nearly indepen- 
dent of Lcxt, down to Pcclct numbers as small as Pe = 4.5. 
In this case, diffusive motion is clearly visible, see movie 
(ring_lowerPE.avi) in the Supporting Material H. This 
implies that the directed two-lane traffic is a very robust 
phenomenon, which is sustainable even in the presence of 
large thermal noise. 

Fig. m indicates that the colloid velocity in the zero- 
density limit is lower for the highest Pcclet number than 
for the intermediate Pcclet number. We attribute this 
behavior to the small but finite rotational Reynolds num- 
ber RCiot — P^^cx)0'^/(4??), which increases with increas- 
ing torque (where p is the fluid density) and reaches 
Rerot = 0.37 for the largest torque; this inertia effect im- 

The movie ring_lowerPE.avi shows a simulation animation of 
2D colloids in a ring channel of median radius R= 2.4(t and width 
D = 2.ba. A constant external torque Lext = 16.7fcsT is applied, a 
factor 9 smaller than for the movie ring_chaiinel.avi. 



plies a small lift force (in analogy with Refs. [Tl[5n]) and a 
reduction of the colloidal density near the walls, in agree- 
ment with the simulated density profiles. 

In order to obtain a net transport, symmetry breaking 
is necessary. The curvature of a ring channel fulfills this 
requirement. We characterize the geometry of the annulus 
by its median radius R and channel width D. We mea- 
sure the average tangential velocity wtan of all colloids in 
the channel for various R and D, keeping the line den- 
sity Nc(t/{2ttR) fixed, where Nc is the number of colloids. 
Results are shown in Fig. [5] for D ~ 1.5(T, D = 2a and 
D = 2.5cr as a function of the median channel radius R. 
We observe a total counter-clockwise tangential motion in 
the annulus that decreases monotonically with increasing 
annular radius for all investigated channel widths. Over 
the accessible range of annular radii, the tangential veloc- 
ity VtioL is well described by a power-law decay, 

Vt.n - R-^. (1) 

Due to the complex interplay of geometric and hydrody- 
namic effects, Vts^niR) does not obey a unique power law. 
For the smallest channel width, D = 1.5, we find 7 ~ 1.7, 
for D = 2.0a an exponent 7 ~ 0.76, and for wider chan- 
nels, D = 2.5a and D = 3.0a, where the colloids can pass 
each other, 7 = 1. This indicates an universal exponent 
for D > 2.5a. 

We propose that the mechanism of net colloid transport 
in a ring channel and its direction can be understood as fol- 
lows. Consider a single colloid spinning counter-clockwise 
about a fixed axis close to one of the confining walls (see 
Fig.[H|). There are three relevant contributions to the force 
acting on the colloid: (i) The surface motion of the colloid 
causes a high velocity gradient in the small gap; due to 
the corresponding shear stress, a force is exerted on the 
colloid parallel to the wall in the direction expected for a 
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Fig. 7: (Color online) Average tangential velocity Vtan of the 
colloids as a function of channel width D for fixed annulus radii 
R = 2.4cr (□, red), R = 3.6a (•, green) and R = 4.8a (A, blue) 
for a constant line density NcO /2nR — 0.8. The torque on each 
colloid is I/oxt.o- 



rolling motion, (ii) A smaller velocity gradient is induced 
in the wide gap, causing a smaller force in the opposite 
direction, (iii) Fluid is pumped into or out from the small 
gap; the resulting pressure difference on both sides of the 
colloid drives it in the direction opposite to the direction 
expected for a rolling motion p^l21) . 

For a single planar wall, these three contributions cancel 
exactly. In case of a linear channel, the contributions (ii) 
and (iii) dominate over (i), causing the two way traffic. 
However, in a ring channel, due to the wall curvature the 
dynamic pressure changes compared to the linear channel. 
Near a convex wall, the flow caused by the spinning colloid 
can avoid the wall more easily than near a concave wall. 
Thus, the pressure gradient near the inner wall is smaller 
than near the outer one. Therefore, the total force acting 
on a colloid near the outer wall is enhanced, and the force 
exerted on a colloid near the inner wall is reduced. This 
explains the difference in the absolute value of the velocity 
near the confining walls that gives rise to the observed net 
transport in the counter-clockwise direction. 

Moreover, the circumference of the outer wall is larger 
than of the inner wall, so that — due to geometric reasons 
— there are less colloids close to the inner wall that could 
diminish the average velocity. The latter effect becomes 
more important at large D and small i?, where the ratio 
of the circumferences of the confining walls increases. 

In a ring channel, colloidal transport can also be 
achieved in very narrow channels, where the colloids can- 
not pass each other, since the forces acting on colloids 
near opposite walls are opposite in direction, but different 
in magnitude. 

When the channel width D is varied for a fixed median 
radius, the tangential velocity shows an interesting, non- 
monotonic behavior, see Fig. [71 For narrow channels, the 
average volocity increases with increasing D and reaches 



a maximum at Z) w 2.1a, where the colloids just cannot 
pass each other. With further increasing channel width, 
the average velocity drops to a minimum at D ^ 2Aa, 
and then rises again. 

In the narrowest channels, with D = 1.5a, we observe a 
low average tangential velocity of the colloids. One reason 
for this is the strong confinement that leads to a reduced 
spinning frequency compared to wider channels. More- 
over, the colloids tend to stay close to the outer wall. This 
configuration allows an unhindered circulation of the fluid 
around the inner confining cylinder in opposite direction 
(compared to the motion of the colloids). 

A maximum in the tangential velocity is reached in 
channels with D 2a for all investigated median radii 
(2.4cr < R < 9.6a). At the maximum, because of the 
additional LJ potential, two colloid just cannot pass each 
other. Since the forces acting on colloids close to the inner 
and outer walls are in opposite directions, but different 
in magnitude, pairs of colloids become wedged together 
into the channel (see snapshot in Fig. [T]) and perceive 
a net force that drives these pairs in counter-clockwise 
direction. As they cannot pass each other, all colloids 
move in the same direction. In this way, the colloids are 
forming plugs that the fluid can hardly pass, the fluid is 
efficiently dragged with these pairs, and a strong back- 
fiow of the fluid that would counteract the motion of 
the colloids is prevented. In the range of D = 2.0a to 
D ~ 2.2a, we observe for R = 2.Aa the formation of col- 
loid triplets as well as pairs and combination of both (see 
movie ring_channel.avi), depending on the initial con- 
ditions, where both types of clusters are stable for roughly 
500 full circulations around the annulus. Since all involved 
interactions are repulsive, the stabilization of the clusters 
is purely hydrodynamics. The highest tangential veloc- 
ity is observed for triplets. For R = 3.6cr and R = A.%a, 
triplets can be observed transitionally, but pairs of colloids 
are the preferred conflguration. 

In wider channels, with D > 2.2a, where the colloids 
can pass each other (see movie ring_channel.avi), col- 
loids at the inner and outer walls move in opposite direc- 
tions, similar to the two-way traffic in the straight chan- 
nels discussed above, and the tangential velocity drops to 
a minimum. However, due to a difference in the absolute 
values of the colloid velocities at the opposite walls, the 
two-way traffic is superposed with an overall translational 
motion and the minimum net velocity is non-zero. With 
further increasing channel width, the tangential velocity 
rises again. This observation is in accordance with the 
case of straight channels, where the average velocity per 
half-channel is higher for broader channels. Moreover, the 
ratio of the circumferences of the two confining cylinders 
increases with increasing channel width, so that there are 
more colloids moving in counter-clockwise direction. Be- 
cause the line density is kept constant, the area density of 
the coloids will decrease with increasing D, implying that 
the density peaks at the walls decrease. As the transport 
is mainly generated near the walls, we expect the veloc- 
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ity to decrease again with further increasing D. However, 
the inner confining cyhnder vanishes for D = 2i? and the 
annulus loses its channel-like character well before. In the 
limit D = 2R, the ring becomes a circular cavity; in this 
case, colloids arc known to move in counter-clockwise di- 
rection tangential to the wall, and experience a small lift 
force HH]. 

Our results in Fig. [7] can be used to estimate typical 
transport velocities of colloids in a ring channel. For the 
experimental system of Ref. [3], $7 = 125 Hz and cr = 
10/im, we find tangential velocities Wtan in the range from 
lOAims"^ to ISO^ims-^ This is more than an order of 
magnitude larger than the transport velocities observed 
for a different channel geometry in Ref. |3]. 

In summary, our simulations show that spinning colloids 
show a rich variety of non-equilibrium states in microchan- 
ncls, from lane formation in straight channels to net trans- 
port in ring channels. In more complex channel networks, 
the fluid flow could be manipulated very easily by con- 
trolling the positions of some colloidal particles by laser 
tweezers, and to use them as valves or to change bound- 
ary conditions. For example, net transport in a straight 
channel could be induced by fixing colloidal particles on 
one side of the channel in regular intervals. Indeed, differ- 
ent kinds of boundary conditions arc another interesting 
possibility to break the symmetry and induce flow, as has 
been pointed out recently for clcctrohydraulic pumping of 
dipolar molecular fluids (like water) in rotating electrical 
fields [22]. The manipulation of dipolar colloids by ro- 
tating electrical fields [5| is another interesting possibility, 
and our results should carry over to such systems to a 
large extent. 

* * * 

Stimulating discussions with C. Bechingcr (Stuttgart) 
arc gratefully acknowledged. 
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